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ABSTRACT 

We present 24 /im images of three protoplanetary disks being photoevaporated around high mass O 
type stars. These objects have "cometary" structure where the dust pulled away from the disk by the 
photoevaporating flow is forced away from the O star by photon pressure on the dust and heating and 
ionization of the gas. Models of the 24 /im and 8 /im brightness profiles agree with this hypothesis. 
These models show that the mass-loss rate needed to sustain such a configuration is in agreement 
with or somewhat less than the theoretical predictions for the photoevaporation process. 
Subject headings: circumstellar matter — stars: protoplanetary disks — stars: formation 



1. INTRODUCTION 



Th e discovery of "proply ds" in Orion fe.g. lO'dell et al.l 
Il99.1t lO'dell fc WeiTll994ft let us image directly proto- 
planetary disks and study their structure. It also led 
to a number of theoretical studies of the photoevapora- 
tion of dis ks by external radiation fields from neighboring 
hot stars I.Tohnstone et alfjflSa IRichling fc YorkfMTflM 
20001: iMatsuvama et all 120031; iHollcnbac h fe Adams! 



2004HThroop feBallvll2005(K Many of these studies pre- 
dict behavior that agrees well with the proplyd proper- 
ties. The proplyds have been observed intensely, for ex- 
amp le through detect ion of associated molecular hydro- 
gen ijChen et alJl998j) and high resolut ion measurements 
in the mid-infrared JSmith et ai]l2005(l . 

Deeper infrared observations can provide additional 
constraints for proplyd models. At 24 /im, the resolution 
of Spitzer is sufficient to identify and perhaps resolve 
such objects at distances up to a few kpc. At longer 
wavelengths, the resolution of Spitzer is too low to re- 
solve these objects, and they will be lost in the confusion 
of the surrounding emission by interstellar material. 

The predicted characteristics of an evaporating disk 
in the ~2 to 50 AU 24 /xm-e mitting zone dif f er sub - 
stantially among the st u dies llJohnstone et al. (1998); 

200j|; IHollenbach fc Adams! J200I : 
(2003). However, the calculations do 



Richling &_Yorl« 



Matsuvama et al 



generally agree that the process occurs on a time scale 
of ~ 1 — 3 x 10 5 years. As a result, one might expect 
to observe a wide range of consequences of disk evap- 
oration. For stars that have dwelt close to an O star 
since their formation, we might find severely truncated 
disks and possibly reduced output in the 24 /im band. 
However, if the stellar velocities within a young stellar 
cluster are random, a star might enter the intense ra- 
diation field of an O star and begin the erosion of its 
disk well after its formation. A system that has just en- 
tered an O star radiation field is likely to be anomalously 
bright at 24 /im because dust grains that would normally 
be shielded from the ambient ultraviolet photons in the 
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optically thick disk are being torn fro m it by the evap- 
oratin g gas and exposed to this field. |Richling fc Yorkel 
( 1998) calculate that the nominal disk emission at 24 
/xm may be boosted more than an order of magnitude 
for some stages in the evaporation process. This paper 
describes three systems that are in this dramatic process 
of early disk evaporation. 

2. OBSERVATIONS 

Our MIPS observations are taken in scan map mode 
and typically cover an area of 0.5-1 square degree. The 
raw data are converted to calibrated imag es with the in- 
strum ent team Data Analysis Tool (DAT) ({Gordon et alJ 
l2?ffl5T) . Point source detection limits are as faint as ~0.4 
mjy, although they can be higher in regions with signifi- 
cant variations in the infrared background. More details 
on the methods for reducing the data (including a de- 
scription of the IRAC data reduction) can be found in 
papers describing the ensemble of infrared sources in the 
three clusters that are the s ubject of this paper: IC 1396 
(jSicilia-Aguilar et alJl2006j) . NGC 2244 (Balog et al. in 
preparation), and NGC 2264 l|Teixeira et alJl2006[) . 

In our full GTO program, we have mapped several re- 
gions including a total of about twenty O stars. Based 
on a hypothesis that an evaporating disk might have a 
cometary appearance at 24 /im, with a photon-pressure- 
driven tail extending radially outward from the O star, 
we inspected these twenty regions. Three clear examples 
were found, illustrated in Figure 1, with their positions 
and integrated flux densities given in Table 1. We found 
no additional sources with tails extending in other direc- 
tions. A simple and very conservative estimate of the 
probability of the three sources being chance alignments 
can be based on the fact that their tails are all within 
10° of radial. This calculation would imply a probabil- 
ity of only about 10~ 4 that the alignments are random. 
This estimate ignores the unique morphology of these 
three sources, so the likelihood is even lower that they 
are randomly occurring phenomena. 

The position of the head of each "comet" coincides 
with a 2MASS and IRAC point source. The sizes of 
the objects are similar to the cometary cl ouds of pho- 
toevap orating globulae found in Carina by iSmith et alJ 
(2004) but the near IR point sources at the positions 
of the "comet" heads indicate that they contain stars. 
No extended structures are observed at wavelengths 
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Fig. 1. — The three evaporating protoplanetary disks together 
with the modeled 24 /im image for the source in NGC 2244. The 
orientation of each image is north up and east toward the left. The 
upper right panel is the model 24 /im image for the tail in NGC 
2244. 



shorter than 24 /im. The only exception is the source 
in NGC2244 where the 8 /im image also shows a simi- 
lar cometary feature. The brightness (Table 1) and the 
infrared color of the near infrared objects in the heads 
indicate that they are stars with masses around 0.5 M 
- 1.0 M Q . These mass estimates probably have large er- 
rors because a significant spread of ages is measured in 
each cluster. 

Table 1. shows the basic parameters of the systems: 
the cluster to which they belong, the position of the head, 
the 2MASS K magnitude of the source in the head, the 
ionizing source and its age, the measured 8 and 24 /im 
flux densities including the tail, the length of the tail, 
the calculated total dust mass and the mass loss rate. 

3. ANALYSIS AND MODELING 

Our analysis is based on the hypothesis that the 
"cometary" tails represent dust entrained in gas evap- 
orating from the working surface of the protoplanetary 
disk and being ejected aw ay from the nearby O-type star. 
iRichling fc Yorkel |l998) show theoretical calculations of 
the early stages of the process that are compatible with 
this possibility. Once the dust has been pulled from 
the disk, the smallest particles (radii <c0.01/im) will be 
stochastically heated and some may be destroyed quickly. 
Small particles (radii < 1000 /jm) will be blown away 
from the O star by photon pressure, or dragged away by 
the gas and will re-emit the absorbed UV energy in the 
thermal infrared. 

The mid-infrared brightnesses of our three objects are 
similar to those of the Orion proplyds, corrected for dis- 



tance ijSmith et al.ll2005|) . Only the system in NGC 2244 
shows slightly extended structure at 8 /im; the other two 
are point-like in the 8 /im images. The fact that all 
the tails are much brighter at 24 /im suggests that small 
grains (size less than 0.1 /im) are the dominant emitters, 
assuming that the grains reach their equilibrium temper- 
atures. To understand the structure of these tails, cuts 
with widths of 6'.'23 were made at 8 and 24/iin along the 
symmetric axis of the tail centered at the O star. The 
background values were estimated using the area above 
and below the tail and subtracted off. The surface bright- 
ness profiles are displayed in Figures El and 0] 

3.1. Basic Model and Parameters 

To model the e mission of the co metary tails we use the 
same model that ISu et, all l)2005t) used for the extended 
Vega debris disk. This model considers the behavior of 
dust grains being blown away from a luminous central 
star by radiation pressure to form an extended outflow 
disk. It should apply generally for outward flows with 
a constant number density radially (n(r) ~ constant). 
To adapt the model to the current situation, we assume 
that the cometary outflow emerges from the disk (comet 
"head" ) position. The remaining parameters are the dust 
properties (composition and grain sizes), the outer ra- 
dius of the tail, and the angle that the line joining the O 
star and the head along the tail makes with respect to 
the plane of the sky (0° means the tail is pointed away 
from the O star exactly on the plane of sky). The emis- 
sion from the tail is computed based on the equilibrium 
dust temperatures assuming the tail is optically thin (no 
stochastic heating is included) . A series of model images 
was computed using different parameters, convolved with 
the Spitzer PSFs to match the observed resolutions, and 
then we extracted radial cuts for comparisons. 

We first explored how the observed profiles of the tails 
vary with different grain sizes and angles to the plane of 
the sky. To simplify the modeling, we assumed the grains 
are astronomical silicates and have a uniform size (radius 
of a) . The modeled 24 /im profiles basically fade rapidly 
as the grain sizes increase. This is because large grains 
have colder temperatures than small ones if they are at 
the same distance from the radiating source; as a result, 
the expected emission at 24 /jm drops significantly, mak- 
ing a shorter observable tail. A similar characteristic is 
seen while varying the angles with respect to the plane 
of the sky, because even though the projected distance 
(on the plane of the sky) remains the same, the actual 
distance between the tail and the O star increases as the 
angle increases, resulting in colder dust temperatures and 
a shorter observable tail. In general, this simple model 
is degenerate between the sizes of grains and the actual 
angles of the tail with respect to the plane of the sky. 
However for the three systems discussed here, models 
based on grains larger than ~1 /im come short of the 
observed tail lengths, and hence cannot fit the data even 
assuming the tails are on the plane of the sky. Without 
other data to constrain the model parameters, hereafter 
we assume the observed tails are exactly on the plane 
of the sky, for simplicity. This assumption is reasonable 
given the selection bias that they would be more difficult 
to recognize at other angles. 

The Tail in NGC 2244: The tail in NGC 2244 is pow- 
ered by the 05 V star HD46150, which we represent with 
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TABLE 1 

The physical and model parameters of the evaporating disks and tails 



Cluster 


Position 


K 


inoizing 


age* 


^8 


F 2 4 


Al 




K 




Aft 




[J2000.0] 


[mag] 


source 


[Myr] 


[mJy] 


[mJy] 


[PC] 




M 




M Q /yr 


NGC 2244 


06:31:54.68 04:56:25.0 


13.53 


HD46150 


4 


1.5±0.2 


77±10 


0.22 


3x10" 


s - 9xl0 -6 


1.4x10 


-9- 4.4 Xl0~ 7 


IC 1396 


21:38:57.09 57:30.46.5 


13.18 


HD206267 


3 


2.3±0.1 


23±5 


0.21 


5x10" 


9 - 2xl0~ 6 


2.4x10" 


" 10 - 7.7 xlO" 8 


NGC 2264 


06:41:01.92 09:52:39.0 


13.65 


S Mon 


4 


2.5±0.1 


78±10 


0.12 


4x10" 


s - lxlO" 5 


3.4x10 


- 9 - 1.1 xl0- e 



*HD46150 IQgura fe Ishida|[T98ll) .HD206267 IHeske. A.^ fc Wcndkcr 19850, S Mon IWalkeJl956l) 

^Dust mass ranges from the observed dust mass using small (~0.01 fim) grains based on the 24 fim profile (lower bound) and total mass integrating 
up to 1 mm size (upper bound). 
^Mass loss rate is computed assuming a velocity of lOkm/s and a gas-to-dust mass ratio of 100. 
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Fig. 2.— The brightness distribution of the tail in NGC 2244 
with diamonds for the data at 24 (im and squares at 8 fim, centered 
at the position of the 05 V star. The contribution of the O star 
has been taken out by PSF subtraction. The 8 ^m profile within 
8" from the O star is contaminated by the residual of the O star. 
The best-fit model is shown with the solid lines. 



a blackbody of 42,000 K and 12 R Q located at a distance 
of 1.5 kpc. As shown in Figure ^ among the three sys- 
tems this tail is the closest to the O star with a projected 
distance of —0.1 pc (therefore the warmest). A tail com- 
posed of a=0.03 fim grains produces a simultaneously 
good fit to the profiles at 8 and 24 /im. Models with 
other grain sizes can also produce good fits to the 24 
fim profile, but the expected emission at 8 fim is either 
much lower than the observed flux (for 0.03 < a < 0.3), 
or much higher for (a < 0.03 ^m). The best fit profiles 
are shown in Figure as solid lines. The tail extends 
from 0.1 to 0.32 pc (Al ~ 0.22 pc) from the O star, with 
a total observed dust mass of — 3xl0 -8 M©. The total 
model emission of the whole system is —74 mJy at 24 
fim and —0.5 mJy at 8 fim. 

The Tail in IC 1396: The distance of IC 1396 is -835 
pc, much closer than NGC 2244. However, the projected 
distance between the disk and the 06 V star HD206267 
is the largest (0.35 pc) with the longest observed angular 
tail length (0.35-0.56 pc, Al — 0.21 pc) among the three 
systems. The tail also appears be clumpy (see Figure GJ, 
although this could also be due to a chance alignment 
with an unrelated source. As in the NGC 2244 object, 
a tail composed of a <0.01 fim grains can provide satis- 
factory fits to the observed 24 fim profile, but it cannot 
reproduce the flux at 8 fim. One possible explanation 
is there exists a population of very small grains (a —10 
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Fig. 3. — The brightness distribution of the tail in IC 1396. 
Symbols used are the same as Figure |21 



A) in the disk that are stochastically heated; therefore, 
the resultant emission spectrum would shift to shorter 
wavelengths producing more 8 fim flux than the emis- 
sion spectrum using the thermal e quilibrium temperature 
(e.g., Figure 8.12 in iKruegell EOQS). The dust mass con- 
tribution from this very small grain population should be 
small. Therefore, the dust mass derived from our simple 
model, — 5xl0 -8 M©, is a reasonable estimate. The to- 
tal model emission at 24 fim is —19 mJy and <0.01 mJy 
at 8 fim. 

The Tail m NGC 2264: The cluster NGC 2264 is lo- 
cated —800 pc away, similar to IC 1396. The star power- 
ing the tail, S Mon, is slightly cooler than the other two 
cases, 07 Ve; therefore a blackbody of 38,000 K and 9 
R© was used. The disk is also farther away from the O 
star, 0.3 pc, with the shortest tail length (0.32-0.44 pc, 
Al — 0.12 pc) among these three systems. In addition, 
the contrast between the head (profile within 6" from the 
head position) and the tail is much larger than the other 
two cases (see Figure 0}. Therefore, similar models used 
to fit the tails in NGC 2244 and IC 1396 cannot obtain a 
good fit to the observed 24 fim profile. An additional un- 
resolved point source at the position of the disk is added 
to fit the profile. Our models are also too faint at 8 fim. 
The behavior at both wavelengths suggests the presence 
of a population of very small grains that are rapidly de- 
stroyed. The models with grains with a <0.01 fim can 
produce good fits to the rest of the 24 fim profile. The 
total observed dust mass is — 4 x 10~ 8 M©, and the total 
emission at 24 fim is —66 mJy and <C0.1 mJy at 8 fim. 
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Fig. 4.— The brightness distribution of the tail in NGC 2264. 
Symbols used are the same as Figure |21 



3.2. Implications 

Although small grains (a ~0.01 /im) are responsi- 
ble for most of the emission at 24 /imi, larger "invisi- 
ble" grains may contribute significantly to the total dust 
mass in the tail. To estimate the mass, we assume that 
the grains have a size distribution similar to the inter- 
stellar medium, N(a) oc a~ 35 , as is also expected for 
grains undergoing collisional cascades in pr otoplanetary 
disks Ipohnanvi I196& iTanaka et al.lll996[L The total 
dust mass in the tail is then M tot oc J ^ira 3 N(a)da oc 



5 da 



,0.5 



Integrating between 0.01 /im and 1 



mm (roughly the largest grain size that can be ejected 
by the radiation pressure of an O star) gives an upper 
limit to the total dust mass in the tail; values are listed 
in Table 1. 

We can also estimate the mass loss rate in the pho- 
toevaporating disks given the velocity of the fl ow in the 
comet ary tails. Following the equation given in lSu et alJ 
(2005), the terminal velocity ranges from ~10 to 100 km 
s _1 for a typical O-type star (luminosity of 10 6 L© and 
mass of 50M Q ) and grains with sizes of 1-0.01 fim, as- 
suming negligible gas drag effects. The presence of gas in 
the tail will result in lower velocities. However, heating 
and ionization acceler ate gas to speeds of 3-20 km s _1 
{Johnstone et aUl998tlBallv et alll99Sfr and these flows 
entrain dust. 

Although the flow velocity is hard to determine accu- 
rately, our estimate for the mass loss rate is robust be- 
cause it depends on the product of the velocity and the 
gas to dust ratio which are roughly inversely correlated 
to each other. For a normal gas to dust ratio (100:1) the 
bulk velocity of the dust and gas mixture is on the order 
of 10 km s _1 . If the gas to dust ratio deviates from the 
standard value (the effect of gas decreases) the bulk ve- 
locity will approach the terminal velocity. Although the 
outflow velocity of material will increase the total calcu- 
lated mass loss will remain approximately the same since 



the lower gas to dust ratio will result in lower total mass, 
hence a lower mass loss rate. 

The X-ray observations of these O stars indicate 
significant stell a r win ds (L x =10 32 ' 4 erg/s for S Mon, 
Simon fc Palm] 12003 L a; =10 32 - 5 erg/s for HD46150, 
Walborn et all l2002t and a brigh t ROSAT source for 
HD206267. iBerghoefer et al.lll996j) . The strong stellar 
winds would speed up the outflow from the disk if the 
disk was very close to the O star (less than 0.01 pc). 
However, such winds should be only a secondary effect 
at the distances of these objects from the stars. 

We can place lower and upper limits on the mass loss 
rate, corresponding to the observed grain sizes and to 
the full size range up to 1000 /im . We take a typi- 
cal gas-to-dust mass ratio of 100 and a velocity of 10 
km/s. The range of mass loss rate is 1.4 x 10~ 10 to 
4.4 x 10~ 8 Moyr- 1 , 2.4 x 10" 11 to 7.7 x 10~ 9 M^r" 1 , 
and 3.4 x 10~ 10 to 1.1 x 10~ 7 Moyr" 1 for the tails in 
NGC 2244, IC 1396 and NGC 2264, respectively. To first 
order, the upper limits of the derived mass loss rates are 
somewhat less than the predictions by the photoevap- 
oration models, - 1 - 6 x 10~ 7 M^yr' 1 at 0.3 - 0.1 
pc from the UV energy source (Richling fc Yorkelll998t 
lYorkd l2004[) . and the lower limits are much less than 
these expectations. 



4. CONCLUSIONS 

We report the detection of three photoevaporating 
disks with cometary tails around high-mass stars in 
Spitzer /24 /im images. Their observed brightnesses and 
morphologies agree well with a model where the dust 
is entrained in the evaporating gas and pulled from the 
disk. Once exposed to the O star radiation field, a struc- 
ture like a comet tail grows, with dust and gas streaming 
away from the O star under the influence of its radiation 
pressure on the dust and heating and ionization of the 
gas. We have roughly estimated the evaporation mass 
loss rate from our observations and find an upper limit 
in reasonable agreement with or somewhat lower than 
theoretical expectations. Further observations of these 
systems should be able to improve our understanding of 
the evaporation process substantially. 
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